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Detection of VER J201 9+407 (G78.2+2.1) and Tycho's Supernova 
Remnant with VERITAS 
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Supernova remnants (SNRs) are widely considered the most likely source of cosmic rays below the knee (10^^ 
eV). Studies of GeV and TeV gamma-ray emission in the vicinity of SNRs, in combination with multi-wavelength 
observations, can trace and constrain the nature of the charged particle population believed to be accelerated 
within SNR shocks. They may also speak to the diffusion and propagation of these energetic particles and to the 
nature of the acceleration mechanisms involved. We report here on the discoveries and interpretation of VHE 
gamma-ray emission from G120.1-I-1.4 (Tycho's SNR) and from the northwest shell of G78. 2-1-2.1 (gamma-ray 
source VER J2019+407, which was discovered as a consequence of the VERITAS Cygnus region survey). 



1. Introduction 



3. Tycho's SNR 



Cosmic rays with energies up to 10^'"' eV are be- 
lieved to be accelerated at strong shocks of supernova 
remnants (SNRs). Current space-borne (Fermi) and 
ground-based (VERITAS, HESS, MAGIC) gamma- 
ray telescopes allow us to study cosmic-ray acceler- 
ators via secondary gamma-ray production, either via 
the inverse Compton process (electrons) or via the de- 
cay of neutral pions (tt" — >■ 77) produced when high- 
energy cosmic-ray ions collide inelastically with the 
ambient medium. Studies of GeV and TeV gamma- 
ray emission in the vicinity of SNRs, in combination 
with multi-wavelength observations, can trace and 
constrain the nature of the charged particle popula- 
tion believed to be accelerated within SNR shocks. 
They may also speak to the diffusion and propagation 
of these energetic particles and to the nature of the 
acceleration mechanisms involved. We discuss here 
recent detections of two TeV gamma-ray sources by 
VERITAS, one identified with Tycho's SNR, the other 
potentially identified with SNR G78. 2+2.1, that have 
the potential to further illuminate this long-standing 
question. 



2. VERITAS 



The VERITAS instrument is an array of four 12-m 
imaging atmospheric Cherenkov telescopes located at 
the Fred Lawrence Whipple Observatory in southern 
Arizona. Each telescope is equipped with a 3.5° field- 
of-view (FOV). The instrument has good (15-25%) en- 
ergy resolution over an energy range of 100 GeV - 30 
TeV, excellent angular resolution (Rgs < 0.1°), and 
achieves a sensitivity of 1% of the Crab Nebula flux 
in approximately 25 hours. For more details on the 
VERITAS instrument and analysis techniques, please 
refer to [Ill,[ll. 



The SNR, G120. 1+1.4, also known as Tycho's 
SNR, is the relic of a Type la supernova[l3| ob- 
served in 1572. Tycho's SNR is young among Galac- 
tic SNRs (438 years). It shows a distinct shell- like 
morphology in radio [5] and strong non-thermal X- 
ray emission concentrated in the SNR rim, with fil- 
aments that have been taken as evidence for electron 
acceleration!! i, [H [11] . 

VERITAS observations performed between 2008 
and 2010 reveal one of the weakest sources yet de- 
tected in TeV gamma rays, with an integral flux 
above 1 TeV of 1.87 ± 0.51stat x lO-i^cm-^s"!, or 
~ 0.9% of the Crab Nebula emission above the same 
energy [l|. The TeV photon spectrum can be described 
by a power law: dN/dE = C(E/3.42TeV)-i', with 
r = 1.95±0.51stat±0.30,ys and C = (1.55±0.43stat± 
0.47sys)xlO-iW-2s-iTeV-i. 

The source is well-fit by a simple symmetric Gaus- 
sian with width fixed at the instrument point-spread 
function (68% containment radius of 0.11°) and is thus 
compatible with a point source. The center of the fit, 
at OOh 25m 27.0s, +64 10' 50" (J2000), is offset by 
0.04° from the center of the remnant, with the sta- 
tistical and systematic uncertainties in this position 
being 0.023° and 0.014° respectively. Figure [T] shows 
that the peak of the gamma-ray emission is slightly 
displaced in the direction of a molecular cloud that 
is seen by integrating over the velocity range between 
^ 68kms~^ and ~ 50kms~^p7| and that may be in- 
teracting with the northeast quadrant of the remnant. 
While it is tempting to interpret this as indicative of 
hadronic emission associated with the cloud itself, not 
only is the effect not statistically significant but the 
association of the cloud with the remnant has been 



recently disputed [23 1. 

Taken together with observations made at other 
wavelengths, however, the case for hadronic acceler- 
ation within the remnant appears strong. Figure [2] 
shows a pair of simple model fits, based on scenar- 
ios where the TeV emission is dominated by leptonic 
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Figure 1: Background-subtracted VERITAS gamma-ray 
count map showing the region around Tychos SNR, 
smoothed with an 0.06 Gaussian kernel. The centroid of 
the VHE gamma-ray emission is indicated with a thick 
black cross 0]. Overlaid on the image are X-ray contours 
(thin black lines) from a Chandra ACIS observation[^ 
and 12CO emission (J=l-0) from the FCRAO Survey 
(gray hues) 3- The VERITAS PSF is shown as a white 
circle. 



(inverse-Com pton ) emission and hadronic (pion de- 
cay) emission[22| respectively, and generated assum- 
ing no influence from the molecular cloud. Particle 
spectra of the form dN/dE = AE-"e(-^/^<^) have 
been used, with the spectral index a fixed to a com- 
mon value for both scenarios. The cutoff energy Ec is 
permitted to differ as is expected for loss-dominated 
distributions 1]. For the leptonic model the derived 
field is ~ 80/iG and the hadronic model has a lower 
limit of ~ 230/iG. These values are well above both 
the usual ~ 3/zG fields of the interstellar medium and 
the values expected from shock-compression of that 
field in moderate-density environments [6^ . As a conse- 
quence, the VERITAS detection of gamma-rays from 
Tycho can be considered additional evidence for mag- 
netic field amplification with the remnant. 

With the models as shown here (fitted only to 
VERITAS, radio, and X-ray data[l|, scenarios where 
the gamma-ray emission is dominated by leptonic or 
hadronic emission are both a reasonable fit, with the 
preference for hadrons resting largely on the fact the 
field of ^ 80/iG required by the leptonic model is no- 
ticeably lower than the values of 200 — 300/iG derived 
from the thinness of the X-ray synchrotron emission 
rimQ. Along these lines of argument, a more re- 
cent paper, incorporating both Fermi-LAT observa- 
tions of Tycho and the VERITAS spectrum shown 
here, claims that the only model consistent with the 
Tycho broad-band spectrum in this case is one in 
which the Fermi and VERITAS gamma-ray spectrum 
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Figure 2: Emission models fitted to the broadband SED 
(radio, non-thermal X-ray, and VHE gamma-ray 
emission) of Tychos SNR. Upper panel and lower panels 
show a lepton-dominated and hadron-dominated model, 
respectively. The long-dashed curve is the 
inverse-Compton (IC) component; the dashed curve the 
pion-decay emission. The solid line at high energies is the 
sum of these components while at low energies it 
corresponds to synchrotron emission. 

is produced by hadrons [20l [2l| . It should be noted, 
however, that regardless of the model used to explain 
the gamma-ray emission, the total relativistic particle 
energy, which represents a significant fraction of the 
total kinetic energy of the SNR, remains dominated 
by hadrons [H. 



4. SNR G78.2+2.1 and VER J201 9+407 

SNR G78.2-f2.1 is a - 1° diameter SNR at ~ 
1.7 kpc distance. At - 7000 years [13, [11 [11 it 
is considerably older than Tycho, and is thought to 
be in an early phase of adiabatic expansion into a 
medium of fairly low density [l^ . Much of the radio 
and X-ray emission lies in distinct northern and south- 
ern features i25(, 30]. Gosachinskij |7| also identifies a 
slowly expanding H i shell immediately surrounding 
the radio shell which Lozinskaya et al. [l^ suggest was 
created by the progenitor stellar wind. A gamma-ray 
pulsar PSR J2021+4026 is located at the center of the 
remnant. Unlike Tycho, G78. 2-1-2. 1's is large enough 
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Figure 3: VERITAS gamma-ray image of SNR G78. 2+2.1 showing the detection of VER J2019+407. The SNR is 
dehneated by CGPS 1420MHz radio contours (white) [ill; the sohd black ellipse co-incident with VER J2019-I-407 is 
the 95% confidence ellipse for IFGL J2020.0-I-4049; the dashed black circle shows the fitted extent of VER J2019-I-407. 
The star symbol shows the location of the central Fermi pulsar PSR J2021-I-4026. The white circle (bottom right 
corner) indicates the 68% containment size of the VERITAS TeV PSF for this analysis. 



to resolve in gamma-rays. 

The VERITAS Cygnus Region Survey^ revealed 
possible emission near G78. 2+2.1, later confirmed by 
~ 20 hours of follow-up observations. In the follow-up 
data alone, a clear signal with 268 ± 34 net counts is 
detected at a location overlapping the northern rim 
of the remnant, significant at the 7.5 standard de- 
viations after accounting for search trials Mod- 
eling of the source by a two-dimensional Gaussian 
convolved with the VERITAS TeV PSF reveals the 
source to be extended, with an intrinsic extension 
of 0.18° ± 0.03°t„t j:|^;^il^^,and a fitted centroid posi- 
tion of R.A. 20''20'"00.0", Decl. +40°49'12" (J2000)). 
However, the apparent source profile (as illustrated in 
Figure[3]) appears both asymmetric and non-Gaussian. 

Figure |3] shows the gamma-ray excess detected by 
VERITAS in relationship to the SNR as seen at other 
wavelengths. The emission seems to follow an arc- 
shaped structure in the radio continuum contours (as 
seen at 1420 MHz). The brightest part of the emis- 
sion is offset by ~ 0.5° from the Fermi LAT pulsar 
PSR J2021+4026 (IFGL J2021.5+4026), but is co- 
located with the Fermi source IFGL J2020.0+4049 
only 9" away. This region of gamma-ray emission does 
coincide with enhanced thermal X-ray emission that 
is suggestive of shocked H i [ll]; curiously, it lies in 
a void of CO emission [l^. The TeV emission also 
lies near a region of bright [S ll] optical line emission 



within the SNR that is identified as shock heated gas 
based on the [S ii]/Ha line ratio 19). 

The relationship of VER J2019+407 and 
IFGL J2020.0+4049 to G78. 2+2.1, and the na- 
ture of the gamma-ray emission being produced, are 
not yet clear. While VER J2019+407 could be a 
PWN related to PSR J2021+4026, the significant 
angular displacement between the two sources, when 
coupled with the total absence of TeV emission at the 
position of PSR J2021+4026, makes this less likely. 
Likewise, while it is possible that VER J2019+407 is 
the PWN of an unknown pulsar in the line-of-sight 
towards SNR G78. 2+2.1, this would ascribe the loca- 
tion of VER J2019+407 near the [S ll] line emission 
and the enhanced thermal X-ray emission to chance 
superposition. A more straightforward explanation 
of the TeV emission and the features observed in the 
X-ray, optical, and radio continuum would be that the 
gamma-ray emission arises from particle acceleration 
(either hadronic or leptonic) within the SNR shock. 
While molecular clouds are frequently invoked as 
potential sites of hadronic gamma-ray emission, the 
dearth of CO in this region means this is not the 
case here; the target material would have to be H 
I, for which the H i shell surrounding the remnant 
could provide a plausible source. In this scenario, the 
pre-shock density estimates of 1 — 10 cm^^ based on 
the ratio of [S ii] A6716 to [S ii] A6731 line fluxes[l| 
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appear consistent with those extrapolated from the 
TeV gamma-ray flux. 

5. Conclusions 

Detection of gamma-ray emission from Tycho's 
SNR and from the vicinity of G78. 2-1-2.1 has added 
to the VERITAS catalog a pair of gamma-ray sources 
with a strong potential to shed light on cosmic- 
ray origin and its relationship to supernova rem- 
nants. Although the relationship between the point- 
like gamma-ray emission from Tycho's SNR and a 
nearby molecular cloud is not clear, the broad-band 
SED from this remnant favors not only a predom- 
inance of accelerated ions within the remnant but 
a hadronic origin to the gamma-ray emission it- 
self. The extended gamma-ray emission from VER 
J2019-I-407, on the other hand, appears to track part 
of G78. 2-1-2. 1's shell, which because of the remnant's 
large size is easily resolved. While the total lack of 
CO at this position makes it clear we are not dealing 
with a molecular cloud scenario, a hadronic origin to 
the gamma-ray emission from VERJ2019-(-407 is not 
ruled out. Preliminary indications are that the den- 
sity of H I, possibly due to interaction with the wall 
of a cavity blown by the progenitor stellar wind, is 
sufficient to produce the observed gamma-ray flux in 
a hadronic emission scenario. 
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